Optical signatures of quantum phase transitions in a light-matter system 



Timothy C. Jarrett, Alexandra Olaya-Castro and Neil F. Johnson 
Centre for Quantum Computation and Department of Physics, 
Clarendon Laboratory, University of Oxford, Oxford, 0X1 3PU, United Kingdom 

(Dated: February 1, 2008) 

Information about quantum phase transitions in conventional condensed matter systems, must 
be sought by probing the matter system itself. By contrast, we show that mixed matter-light sys- 
tems offer a distinct advantage in that the photon field carries clear signatures of the associated 
quantum critical phenomena. Having derived an accurate, size-consistent Hamiltonian for the pho- 
tonic field in the well-known Dicke model, we predict striking behavior of the optical squeezing and 
photon statistics near the phase transition. The corresponding dynamics resemble those of a degen- 
erate parametric amplifier. Our findings boost the motivation for exploring exotic quantum phase 
transition phenomena in atom-cavity, nanostructure-cavity, and nanostructure-photonic-band-gap 
systems. 
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There are several theoretical models which are cur- 
rently attracting attention, based on the possible insights 
that they offer into the nature of Quantum Phase Transi- 
tions (QPTs). One of these is the Dicke model which was 
originally developed in quantum optics, together with its 
recent generalizations Q, I^Q. In practice, such exotic 
quantum phenomena can only be studied experimentally 
if the system's many-body quantum state can be probed 
in some way T] . Unfortunately in condensed matter sys- 
tems, such probing is typically indirect since one cannot 
'see' many-body quantum spin states. 

Here we predict that in light-matter systems approx- 
imating to the Dicke model ~ such as atom-cavity and 
nanostructure-cavity systems 0, S B ~ the statisti- 
cal properties of the photon field offer direct and strik- 
ing signatures of the quantum critical phenomena un- 
derlying a QPT. Our results are based on an accu- 
rate, size-consistent calculation of the statistical prop- 
erties of the photon field in the Dicke model, and 
help motivate the exploration of such exotic quantum 
phenomena in atom-cavity, nanostructure-cavity, and 
nanostructure-photonic-band-gap systems [H Q, y, |j| . In 
addition to the results themselves, our theoretical ap- 
proach presents a number of distinct features over previ- 
ous works (1) we avoid using canonical pertur- 
bation schemes and projection methods, which can suf- 
fer from inconsistent size-dependencies as one approaches 
the thermodynamic limit ja]- Instead we adopt a simi- 
lar renormalization-like scheme to Ref. but take the 
opposite viewpoint by renormalizing the dynamics of the 
photon field as opposed to the matter system. (2) Our 
approach shows the direct connection between the Dicke 
model and a degenerate parametric optical amplifier. (3) 
In addition to explicitly reproducing the correct scaling 
near the critical point, we are able to show that striking 
signatures arise in a number of key statistical properties 
associated with the photon field. (4) We are able to show 
the optical manifestation of a quasi-integrable to quan- 
tum chaotic transition near the QPT. 

The Dicke model describes the interaction between a 



single-modephoton field and N non-interacting two-level 
systems flQ: 

H = a)a + eJ:, + ^^{a + a))J,^, (1) 
2y N 

where J, = ^ E^i < and = +'^r) are 

the collective angular momentum operators, and the op- 
erators a,a^ and crfjCr^ correspond to the photon field 
and two-level atom i respectively. This model exhibits 
a phase transition at both zero and finite temperature 
Qil^ ■ Employing the cumulant expansion method 0, |^ , 
we here choose to eliminate the degrees of freedom in 
the matter subsystem and hence derive a size-consistent 
effective Hamiltonian for the photon field. Consider 
H = Ho + Hi where 

Ho = Hb = eJz 

2A 

Hi = Ha + Hab = a^a H — (a^ + a) J^, 

with the Hamiltonian Ha denoting subsystem a (i.e. pho- 
ton field), Hb denoting subsystem b (i.e. matter) and Hab 
denoting their interaction. The size-consistent ^ effec- 
tive Hamiltonian in subsystem a is 

The index c denotes cumulant averaging f3| and 
the thermal average is carried out with respect 
to the matter degrees of freedom, i.e. {A)b = 
Trfc (e~f^^-^'A) /Tn {e'f^"-^'). The LiouviUian superoper- 
ator Lq is defined by LoA — [Hb,A]. The cumulants can 
be expanded in a series: 

^^-PiH,+L.) .I)c^f2 ^-Tm + LoYY- 
v—1 

We calculate the first two cumulants exactly: 

{Hi + LoY - a^a 
{{Hi + Lofr = XHa + a^f. 



2 



When we calculate higher-order cumulants, and retain 
terms which are linear and quadratic in a and a\ we 
can derive a general expression for every even and odd 
cumulant: 

{{Hj + Laf'^+^Y = A2e(2"+i)-2(a + a^f tanh (^^ j . 

We find that only second-order terms are required to de- 
scribe the salient features of the Dicke model at zero- 
temperature, hence we will exclude higher-order terms 
from the discussion. The effective Hamiltonian is there- 
fore 



a) a 



A2 



■ tanh 



(a + a)y 



In the low-temperature limit [3 



this reduces to 



H 



a(T=0) 



(2) 



with bj = 1 — ^'^d 7 = — Hence the optical prop- 
erties of the Dicke model can be mapped onto a degener- 
ate parametric process in which a classical field interacts 
with a non-linear medium. Equation (jSJ belongs to a 
class of squeezing Hamiltonians with SU(1,1) symmetry 
which have been shown to exhibit a ground-state phase 
transition 
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FIG. 1: Ground state energy as a function of the light- matter 
coupling A, as calculated with the size-consistent Hamiltonian 
of Eq. (|5J (solid line) and the full Dicke model Hamiltonian 
(circles). 

We now show how our effective Hamiltonian not only 
reproduces the important features of the Dicke model, 
but also provides interesting optical signatures of its un- 
derlying quantum critical behavior: 

(i) Ground state energy, quantum critical point and scal- 
ing. Figure n shows the excellent agreement between the 
exact numerical ground-state energy of the Dicke-model, 
and that obtained from Eq. 0. In order to describe 



the underlying critical behavior, we perform a Bogoli- 
ubov transformation from the operators a,o^ to squeezed 
6-bosons, such that — ua + va^ and |up — \v\'^ ~ I: 



6 = 



vr 



a' 



aa 



and substitute for a,a^ in Eq. The additive constant 
has no effect on the energy gap between the ground and 
first excited states. Using [6, b''] = 1 yields 



H 



eff _ 
a(T=0) " 



1 -a2 



(7a 



-7) {b^' 



(3) 



We choose a such that the coefhcient of {b^^ + b^) is zero 
and apply the so-called resonance condition e — 1. We 
find 



1 - 2A2 ± VI - 4A2 
2A^ 

and the effective Hamiltonian (having arbitrarily chosen 
one of the two equivalent solutions) becomes 



(4) 



This effective Hamiltonian represents a simple harmonic 
oscillator in the 6-bosons and hence may be diagonalized 
by the num ber states of the 6-boson operators, i.e. _E„ = 
-V1-4A2 (n-^i). The energy gap between the ground 
and first excited states is 



AE = Vr^^4A2 

and hence we correctly deduce the quantum critical point 
as A = Ac = 0.5. The energy gap is proportional to 
\/Xc ~ X as A — > Ac, which is also in excellent agreement 
with previous results 0. The sub-radiant phase (i.e. 
A < Ac) is well described by the above Bogoliubov trans- 
formation - the effective Hamiltonian maps to a simple 
harmonic oscillator. By contrast, in the super-radiant 
phase (i.e. A > Ac) the effective Hamiltonian in Eq. Q 
resembles an inverted oscillator, as we discuss later. This 
finding highlights the inherent instability associated with 
the phase transition. 

(ii) Sub-radiant phase. We obtain analytical expressions 
for the photon field occupation number TV, the Mandel Q- 
parameter (?[ and the optical squeezing [3 , by exploiting 
the SU{1, 1) symmetry of the effective Hamiltonian (Eq. 
@) 0]. The Lie algebra oi SU (1,1) is generated by 
introducing the operators 0| 



^ 2 



K_ = —aa, 
2 ' 



with [Ko,K±] ^ ±K±,[K^,K+] = 2Ko. Equation © 
becomes 

=2c^ifo + ^7(^+ + ^-) 
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with 7 redefined as 
coherent states 



We introduce the S'?7(l,l) 



\^,k) =exp{zK+- z*K^}\0,k), 

where z = -{9/2)e-"l' and ^ = - tanh(6'/2)e-*'^; 9 and 
(j) are group parameters with ranges (— oo, oo) and [0, 27r] 
respectively. The state |0, k) with fc = i is the vacuum 
squeezed state while fc = f corresponds to the squeezed 
one-photon state. Now N, Q and the squeezing can be 
calculated in terms of expectation values of the three 
generators. For example 

{^,k\Ko\^,k) = /ccoshe, 
{^,k\K±\^,k) = -ksinhOe'^"''. 

The equations of motion for 9 and </> are 



— (fc sinh( 
(fc sinh 6) 



~d9' 



where Jif — ^Ha{T=o)\^' hence obtain 



1 

2' 



J/f = 2uj k cosh 9 — "fk sinh 9 cos ( 

9 = — 7sin0, 

(j) = 2uj — J coth 9 cos (j) . 



The stationary points must satisfy 9 — (j) = 0, yielding 
(j) — nTT with n = 0, 1, 2, . . .. If n is odd, we find 




FIG. 2: Analytical (circles) and numerical (solid line) calcu- 
lation of the photon occupation number plotted as a function 
of A. The inset shows the saturation of A'^ for large values of 
A. 



Fig. and the ground state is a minimum-uncertainty 
squeezed state, i.e. (AXi)2(AX2)2 = j^. The Mandel 
Q-parameter in the sub-radiant phase is given by 

^ (gtagta) - (ata)2 _ 4(j^g) - 4(j^o)^ 
^ (ata) 2(i^o>-i 

fc((l + 2fc)cosh(26l)-f 2A:-l)-4Pcosh^6l ^ 
~ 2fccosh6' - i 

which, for fc = i, is plotted in Fig. El^c). 



9 — arccoth 



1-2A^ 



Hence no real solution exists above the critical point 
Ac = 0.5. As hinted at earlier with the suggestion of 
an inverted oscillator, this is because the energy Jff be- 
comes unbounded above Ac. In the sub-radiant phase 
(i.e. A < Ac) the photon number N = k\a^ a\^, k) can 
be written as 

N ^ 2{Ko) -^ = 2kcosh9-^ . 

Figure|21shows that this expression accurately reproduces 
the numerical results, implying that the field is effec- 
tively in a squeezed state in the sub-radiant phase. We 
now turn to the optical squeezing itself, which is defined 
in terms of the quadrature operators Xi = ^{a + a') 
and X2 = ^(a - a^) with [Xi,X2] = i/2. This yields 
(AXi)2(AX2)2 > ^ where (AX,)^ = (Xf) - {X,)\ 



Hence squeezing exists if {/S.XiY 



< 



We obtain 



(AXi 



k (cosh 9 =F cos sinh 9) . 



For k = J and (p = nn with n = 0, 1, 2, . . ., we find 
that (AX2)^ < J throughout the sub-radiant phase (see 
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FIG. 3: Analytical (circles) and numerical (solid line) calcula- 
tions of (a) variance of Xi, (b) variance of X2 and (c) Mandel 
Q-parameter vs. A. 



(Hi) Super-radiant phase. We introduce the canonical 
position and momentum operators 



x= ^={a + a^); p^^J'^{a~a^) (5) 



2uj 



4 



where \x^p\ = i. Equation (0) becomes 
with 



trer r c . ^ 2-2 



(6) 



m 



1 _ 22 



-2 2 / 1 , ^7 
w = w 1 H 



1 - 



27 



where we have again ignored the additive ^. The coeffi- 




FIG. 4: CoefRcients of (solid line) and (dashed line) in 
Equation @ vs. A. 

cients of and are plotted against A in Fig. 21 These 
plots show that in the sub-radiant phase (i.e. A < 0.5) 
the system is equivalent to a harmonic oscillator and as 
such is square-integrable - hence it can be diagonalized 
using basis states of the number operator a^a. In the 
super-radiant phase (i.e. A > 0.5) the system becomes 
an inverted harmonic oscillator - first in momentum, and 
then in position. The properties of an inverted potential 
harmonic oscillator are discussed elsewhere It is in- 
teresting to note that the energy becomes unbounded 
and the energy eigenspectrum becomes continuous - in 
addition, the inverted potential harmonic oscillator has 
been used as a model of instability in relation to quantum 
chaos . This confirms the claim of Ref. ^3 that the 
system undergoes a transition from quasi-integrable to 
quantum chaotic behavior at the quantum critical point. 



In the super-radiant phase, the system cannot be de- 
scribed in a simple analytic form - except in the limit 
A — > oo. However the quantities in which we are inter- 
ested can still be calculated using numerical simulations. 
These numerical simulations suggest that - with the ex- 
ception of the energy - each quantity of interest rapidly 
converges to a single value as A is increased. Remark- 
ably, we find that only a small number of photons are 
required to obtain extremely accurate results. Further- 
more, the scaling remains intact as we increase the sys- 
tem size. Figure |21 shows that the occupation number 
N is initially zero, and remains small but finite as we 
increase A toward Ac. At this critical point, the occu- 
pation becomes macroscopic, i.e. the system ceases to 
be sub-radiant and enters the super-radiant phase. The 
behaviour seen in Fig|21 compares extremely well with 
the properties predicted for the output photon flux in 
Ref. Indeed Fig. |21 (inset) suggests that as A in- 

creases further, the occupation converges toward a single 
value. Figures Ola) and (b) indicate that the radiation 
field survives in an ideal squeezed state up to the critical 
point Ac = 0.5, but that this squeezed state then breaks 
down in the super-radiant phase. At Ac, a discontinu- 
ity appears in (AX2)^ which is associated with a sudden 
change in the photon statistics - a finding confirmed by 
the Mandel Q-parameter (FigOIc)). In the sub-radiant 
phase (i.e. A < Ac) we have Q > indicating that the 
photon statistics are super-Poissonian. At Ac, Q diverges. 
Above Ac, Q < which indicates that the statistics are 
sub-Poissonian. As A is increased further, Q —^ —1 and 
hence Q takes on its most negative possible value, i.e. 
the multi-photon state becomes a Fock state 0- 

To summarize, we have analyzed a known Quantum 
Phase Transition from the entirely new perspective of the 
accompanying photon field. The predicted optical signa- 
tures should be readily measurable in either atom-cavity 
or nanostructure-cavity systems 0, Qi 0| using existing 
optical techniques 0- 
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